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Abstract 

li F : H ^ H is a map in a Hilbert space H, F £ Cf^^, and there exists 
y, such that F{y) = 0, F'{y) / 0, then equation F{u) = can be solved by a 
DSM (dynamical systems method). This method yields also a convergent iterative 
method for finding y, and this method converges at the rate of a geometric series. 
It is not assumed that y is the only solution to F{u) = 0. Stable approximation to 
a solution of the equation F{u) = / is constructed by a DSM when / is unknown 
but fs is known, where — /|| < 6. 

1 Introduction 

In this paper a method for solving fairly general class of nonlinear operator equations 
F{u) = in a Hilbert space is proposed, its convergence is proved, and an iterative 
method for solving the above equation is constructed. Convergence of the iterative 
method is also proved. These results are based on the following assumptions: a) the 
above equation has a solution y, possibly non-unique, b) F e C^i^^y and c) F'{y) ^ 0. 
No restrictions on the rate of growth of nonlinearity are made. The literature on methods 
for solving nonlinear equations is large (see, e.g., [T] and references therein). Most of the 
results obtained so far are based on Newton-type methods and their modifications. There 
is also a well-developed theory for equations with monotone operators (0)- The method 
used in this paper is a version of the dynamical systems method which is studied in Pj. 
The idea of this method is described briefly below. 

Let F : H ^ H he a. map in a Hilbert space. Assume that equation 

F(u) = (1) 
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has a solution possibly non-unique, and 

F'{y) + 0, (2) 

where F' is the Frechet derivative of F . This ssumption means that F\y) is not equal 
to zero identically on H. Thus, it is a weak assumption. Assume that F G C^^^^ i.e., 

sup \F^^{u)\ < Mj{R) < i < 2, (3) 

ueB{uo,R) 

where uq E H is a. given element, R > 0, and no restrictions on the growth of Mj{R) 
as R grows are made. This means that the nonlinearity F can grow arbitrarily fast as 
II M — Moll grows. It is known that under these assumptions equation (Q) may have no 
solutions. Thus, we have assumed that a solution y to ^ exists. 

We do not assume that F'{u) has a bounded inverse operator. Therefore the standard 
Newton-type methods are not applicable. Dynamical system method (DSM) consists of 
finding an operator $ such that the problem 

u = ^{t,u), u{0) = uo (4) 

has a unique global solution u{t), there exists u{oo), and F{u{oo)) = 0. To ensure the 
unique local solvability of (0)) we assume that 

\mt,u) - $(t,f)|| < L{R)\\u - v\\ ^u,v e B{uo,R). 

Then the global existence of the unique local solution holds if sup^ ||^(^)|| < oo. 
The results of this paper can be summarized in two theorems. Let us denote 

A:= F'{u{t)), T := A*A, Ta := T + aF, A:=F'{y), f = A*A. (5) 

Assume that a(t) is a positive monotonically decaying function, 

a{t) > 0, lim a{t) = 0, — < J. (6) 

Theorem 1. If a solution y to equation (Q) exists, possibly is non-unique, and assump- 
tions and Q hold, then y = u{oo), where u{t) solves the following DSM problem: 

ii = -T-yA*F(n(t)) + a{t){u{t) - z)], ti(0) = «o, (7) 
and where z and uq are suitably chosen. 

Theorem 2. Under the assumptions o/ Theorem ^ the iterative process 

Un+l = Un- hnT~^^[A*{Un)F{Un) + UniUn- Z)], Uq = Uq, (8) 

where hn > and > are suitably chosen, generates the sequence Un converging to y. 
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Remark 1. The suitable choices of Un and are discussed in the j^roo/ 0/ Theorem |21 

Remark 2. Essentially, Theorem says that any solvable operator equation with Cf^^ 
operator, satisfying only a weak assumption can be solved by a DSM. Condition 
means that the range of the linear operator F'{y) contains at least one non-zero element. 
It allows F'{y) to have infinite- dimensional null-space. 

In Section El we prove Theorem^ and Theorem |21 In their proofs we use the following 
lemmas. 

Lemma 1. Assume that g{t) > is a C^([0,oo)] function satisfying the inequality 

g<-l{t)g + a{t)g^ + m. ^ > 0, (9) 

where 7, a and (3 are nonnegative continuous functions. Assume that there exists n G 
C^([0,oo)), fi> 0, \imt-,oo fJ'it) = 00, such that 

(10) 

Then any solution to exists on [0, 00) and 

0<git)<-^y te[0,oo). (11) 

Lemma H is proved in [21 pp. 66- 70]. 

Lemma 2. Let gn+i < jgn + pgn' (7o > 0, < 7 < 1, p > 0. If m < where 

7 < g < 1, then lim„^oo = 0, and gn < 5'o«?"- 

Proof of Lemma\^ Estimate gi < 7m + pm'^ < qm holds if m < 7 < g < 1. 
Assume that gn < g^q^- Then 

Qn+i < igoq"" +p{goq"f = Qaq" ij + pgoq") < 9oq"^\ 
because 7 + pgoq"^ < 7 + pgoq < q- Lemma 2 is proved. □ 



2 Proofs 

Proof of Theorem^ If F'{y) := A is linear and A ^ 0, then there exists vi 7^ 0, fi = Tv. 
By the linearity of T, every element cvi belongs to the range of T for any constant c, 
because cvi = T{cv). Therefore there exists a z such that y — z = Tv, where > 
can be chosen arbitrarily small. How small \\v\\ should be chosen will become clear later. 
Let u{t) — y := w(t), \\w(t)\\ := g{t). Write equation as 

w = -T.I [A*{F{u) - F{y)) + a{t)w + a{t){y - z)], (12) 
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and use the formula F{u) - F{y) = Aw + K, where \\K\\ < Then 

w = -w- T;^1^ A*K - ait)T-^l^fv. (13) 
Multiply this equation by in if and use the estimate ||Tg7^yl*|| < a > 0, to get 



gg < I M2g' / , _ , , \ 

9./^ 2 V y II \^ a(t) a{t) a{t)J II 



WW. 



If a > then \\f-'^f\\ < 1, 

a||r-i < 1, a\\{T-' - f-')f\\ = a\\T-\A*A - A*A)f-'f\\ < 2M,M,g. 
Collecting the above estimates and choosing II^H so that 2MiM2||t'|| < |, we obtain 

g < -- + ^== + a{t)\\v\\, Co:=— . (14) 
2 ^a{t) 4 

Apply Lemma lU to Jll). Here 7 = i a = -a=, p = a{t)\\v\\. Let ^lit) = A = 

const > 0. Condition i) of Lemma IT] holds if — < — ^= ( ^ — ^ — ). This inequality 
holds if A > 8co, see the last assumption ©• Condition n«) holds if g(0)—^= < 1. 

A/a(0) 

This inequality holds (for any initial value g{0) = \\uo — y\\) if a(0) is sufficiently large. 
Condition ii) holds if y'a(t)||i;|| < where we have used the last assumption © 
again. This inequality holds if (*) 8Aiya(0)||i;|| < 1. (Recall that a(0) > a{t) due to 
monotonicity of a{t).) Inequality (*) holds if \\v\\ is sufficiently small. Thus, if is 

sufficiently small, then, by Lemma [TJ we get g{t) < ^/^.^ so \\u{t) — y\\ < g as 

t ^ 00. 

Theorem ^ is proved. □ 

Proof of Theorem\^ Let Wn '■= Un — y, g-n '■= \\wn\\- As in the proof of Theorem 1, we 
assume 2MiM2||f || < | and rewrite (8) as 

Wn+i =Wn-K [A*{un){F{un) - F{y)) + a„«;„ + an{y - z)], Wo = \\uo - y\\. 

Using the Taylor formula 

F{Un) - F{y) = A{Un)Wn + K{Wn), \\K\\ < 

the estimate ||T^^y4*(w„) || < 2^^, and the formula y — z = Tv, we get 

Wn+l = (1 - K)Wn - hnT„^A*{Un)K{Wn) - KanT'^fv. (15) 
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Taking into account that \\T^ ^T\\ < 1, and a\\T^ ^|| < 1 if a > 0, we obtain 
and 

IK^.:.^ - T-:)f\\ = iiT-;(f.„ - T.jf-;f II < ^-^^i^ ■.= ^ 



Let Co := Then we obtain from p5|) the following inequality: 

gn+l < (1 - K)gn + + Cihn\\v\\gn + hnan\\v\ 



We have assumed in the proof of Theorem 1 that ci\\v II < |. Thus 



9' 



Choose an = IQclgl Then ^ = i, and 



" 2 I 7 II I 

n+l < (1 - -^)9n + —^9n + KcLnM 



9n+i<{l-'-^)9n + l'ocQK\\v\\gl, 9o = \\uQ-y\\<R, (16) 



where i? > is defined in Take ^„ = h G (0,1) and choose := m such that 

) and q + h > 1. Then Lemma |2 

Un — ?/|| < (?o?" — > as n — > cx). 



m < x6c-ofe|H| ' ^^^^^ ^ ^ (0' 1) ^^'^ q + h > 1. Then Lemma El implies 



Theorem |21 is proved. 

□ 



3 Stability of the solution 

Assume that F{y) = /, where the exact data / are not known but the noisy data fs are 
given, 11/5 — /II < S. Then the DSM yields a stable approximation of the solution y if 
the stopping time ts is properly chosen. The DSM is similar to (7): 

its = -T;^l^[A*{Fiusit)) - fs) + ait)iusit) - z)], us{0) = Wq, (17) 

Let ws := us{t) — y, gsif) := \ \ws\\- As in the proof of Theorem 1 we derive the inequality 
similar to (14): 

95 , Cog] , 5 _ M2 
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and apply Lemma 1. The only difference is in checking condition ii) of Lemma 1. This 
condition now takes the form: 



a(t)\\v\\ + — ^= < -^^^-f^. 

This condition can only be satisfied for t G [Ojt^], where ts < oo. The stopping time ts 
can be determined, for example, from the equation = |, provided that v is chosen 

sufiiciently small, so that 8Aa/ a(0)| |v| | < |, and A > Scq as in the proof of Theorem 1. 
Then, by Lemma 1, we have gs{ts) < ^^^^ — > as 5 — > 0. Let us formulate the result. 

Theorem 3. Let us := us{ts), where us{t) solves problem (17) andts is chosen as above. 
Then lim5_»o | |uj — y| | = 0. 

References 

[1] Ortega, J., Rheinboldt, W., Iterative solution of nonlinear equations in several 
variables, SIAM, Philadelphia, 2000. 

[2] Ramm, A. G. , Inverse problems. Springer, New York, 2005. 

[3] Zeidlcr, E., Nonlinear functional analysis and its applications, II-B, Springer, New 
York, 1985. 



6 



